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ABSTRACT 

GD 356 is unique among magnetic white dwarfs because it shows Zeeman-spUt 
Bahner hues in pure emission. The hues originate from a region of nearly uniform field 
strength {5B/B w 0.1) that covers 10 per cent of the stellar surface in which there 
is a temperature inversion. The energy source that heats the photosphere remains a 
mystery but it is likely to be associated with the presence of a companion. Based on 
current models we use archival Spitzer IRAC observations to place a new and stringent 
upper limit of 12 Mj for the mass of such a companion. In the light of this result and 
the recent discovery of a 115 min photometric period for GD356, we exclude previous 
models that invoke accretion and revisit the unipolar inductor model that has been 
proposed for this system. In this model a highly conducting planet with a metallic 
core orbits the magnetic white dwarf and, as it cuts through field lines, a current 
is set flowing between the two bodies. This current dissipates in the photosphere of 
the white dwarf and causes a temperature inversion. Such a planet is unlikely to have 
survived the RGB / AGB phases of evolution so we argue that it may have formed from 
the circumstellar disc of a disrupted He or CO core during a rare merger of two white 
dwarfs. GD356 would then be a white dwarf counterpart of the millisecond binary 
pulsar PSR 1257+12 which is known to host a planetary system. 

Key words: white dwarfs ~ planetary systems 



1 INTRODUCTION the polarization spectra. Barring the possibility of a rapid 

rotator, the observed lack of variability of the spectrum over 

GD 356 was first noted to be peculiar among magnetic white periods of several hours to days must indicate that the spin 

dwarfs through the discovery of strong Zeeman split Balmer ^nd dipole axes are nearly aligned. Likewise, the absence of 

lines in pure emission jGreenstein fc McCarthv| | l985| ). The ^gpo^s of significant changes in the occasional spectra taken 

number of known magnetic white dwarfs has since increased ^^gj. ^ 35 yr period (G. Schmidt and J. Liebert, private com- 

to some 220 through surveys such as SDSS but GD 356 re- munication) suggests that the photospheric region that gives 

mains unique in exhibiting these properties. ^-^^^ to the emission lines has a stable structure on this time 

A detailed stud y of the atmosphere of GD 356 by scale. 
iFerrario et al.l (Il997l ) established that the line and contin- 
uum spectra can be modelled with a white dwarf with an The fields of the high field magnetic white dwarfs 
effective temperature T^s = 7 500K and an assumed grav- (HFMWDs), of 1 < B/MG < 1000, such as GD356 are 
ity \ogiQ{g/cms~^) = 8 with a centred dipole field structure generally believed to be of fossilized origin rather th an dy- 
and polar field strength Bp = 13 MG. The emission lines namo generated (jWickramasinghe fc Ferrarid I2OO0I ). They 
arise from a ring like spherical sector or strip around the are most likely generated in a co mmon envelope th at left a 
magnetic pole covering lOper cent of the stellar surface. The very close binary or merged core (|Tout et al.ll2008l ). In this 
lines are a result of a temperature inversion that begins deep case all single white dwarfs with high magnetic fields have 
within the photosphere at optical depths of 1.0. The remain- evolved from binary stars that merged during common enve- 
der of the photosphere produces absorption lines broadened lope evolution or shortly afterwards. Though very uncertain, 
by the underlying dipolar field. These lines are masked by estimates of their number densities indicate that about three 
the emission lines in the flux spectra but are discernible in times as many systems that enter a common envelope phase 
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of evolution end up merging as form cataclysmic variables. 
The long ohmic decay time scales, 8 — 12Gyr l|Cummind 
of the dipolar component in magnetic white dwarfs 
and the lack of observational evidence for differences in the 
mean dipolar field strength along the white dwarf cooling se- 
quence supports this hypothesis. Cool white dwarfs develop 
convective envelopes that could potentially lead to genera- 
tion of magnetic fields by a contemporary dynamo but there 
is no evidence for an increase in the incidence of magnetism 
among the cooler and more convective stars in the well stud- 
ied high-field group of white dwarfs. It is of course possible 
that all cool white dwarfs with outer convective envelopes 
have dynamo generated fields that are below the current 
observational limit of detectability fro m Zeeman polarime- 
try of B ^ 10^ G (| Jordan et al.l 120071 '). However, attempts 
to detect X-ray emission from coronae that may be gener- 
ated through magnetic activity in such stars have so far led 
to upper limits that are well below theoretical predictions 
iMusielak et al.ll2003l ). 

Since these early investig ations there have been s ome 
new observations of this star. ' Brinkworth et aD (|2004l ) re- 
ported the detection of low-amplitude (±0.2 per cent) near 
sinusoidal variability in the V band with a period of 115 min. 
They attribute this to the rotation period of the star. They 
presented a model in which a dark spot covers 10 per cent 
of the stellar surface and is viewed nearly face on or edge 
on as the star rotates. They speculated that the temper- 
ature inversion required to explain the spectroscopic data 
must occur in this region. While the cause of the observed 
photometric variability may well be related to a temperature 
differential between the line emission region and the remain- 
der of the star, the idea of a dark magnetic spot presumably 
with an enhanced field caused by magnetic activity as in 
the Sun and stable over some 25 yr is less attractive. The 
modelling shows no evidence for an enhanced field strength 
in the emission line region. Rather it is simply a specially 
heated region of the star with an otherwise approximately 
dipolar field structure. It is thus more likely that GD 356 has 
a fossil field like other white dwarfs of similar field strength. 

The source of energy that powers the emission line re- 
gion in GD 356 remains a mystery. The observed luminosity 
in the Balmer lines is 2 x lO^^ergs s~^, much larger than 
the stringent upper limit of 6 x lO^^ergs s~^ that has re- 
cently been placed on the X-ray luminosity of GD 3 56. So the 
source of heating is not an X-ray corona (|Weisskopf et al.l 
I2OO7I ). Given the implausibility of a single star interpreta- 
tion, the most likely possibility is that the energy is ex- 
tracted in some way from a companion. The lack of evidence 
for an accretion disc in the line spectrum or of emission from 
accretion shocks appears to preclude intermediate polar type 
and AM Her type models. 

One of the more intriguing models that has been pro- 
posed for GD 356 assumes that it has a companion of plan- 
etary mass with a conducting compo sition in a close orbit 
jLi. Ferrario fc Wickramasinghdl 19981 ). Such a planet would 
act as a unipolar inductor generating an electrical current 
that fiows between the two stars analog ous to the model pro- 
posed for the Jupiter-Io system (Gol dreich fc Lvnden-BeTl 
1 19691 ). Ohmic dissipation in the atmosphere of the white 
dwarf leads to a temperature inversion in a ring surround- 
ing the magnetic pole as the planet orbits the white dwarf. 
This model, as it presently stands, hinges rather heavily on 



Table 1. Mid- Infrared Fluxes for GD356 



Wavelength 


Model Flux 


Measured Flux 


X/fim 






3.6 


512 


523 ± 26 


4.5 


339 


347 ± 17 


5.7 


218 


222 ± 12 


7.9 


121 


138 ±9 



There is a 2. Oct excess measured at 7.9 ^m. 

the possibility of the survival of such a planet through the 
RGB and AGB phases of evolution. 

In this paper we present an analysis of archival 
Spitzer/IRAC observations of GD 356 that rule out a com- 
panion of mass as greater than 12 Mj , according to substellar 
cooling models. We then discuss possible accretion models 
and the unipolar inductor model in the light of the new mass 
limits and the recently discovered photometric period for 
this star. We argue that the ensemble of presently available 
observations is inconsistent with accretion from a gaseous 
planetary companion being the source of the anomalous line 
emission seen in GD 356. We then show that rocky planets 
are unlikely to survive the AGB phase of evolution and be 
dragged into a close orbit as is required in the unipolar in- 
ductor model. We propose that, for this model to be viable, 
the planet must have formed from material in a disc that re- 
sulted from the merging of two white dwarfs and argue that 
GD 356 may be the result of such a rare event akin to the 
event that resulted in the format ion of planets around th e 
millisecond pulsar PSR 1257-1-12 l|Wolszczan fc Fraillirgg^ ). 



2 THE OBSERVATIONS, LIMITS ON 

COMPANION MASS AND IMPLICATIONS 

We combine multi-wavelength photometry from several 
sources in order to constrain the spectral energy dis- 
tribution of GD 356, particularly the photospheric emis- 
sion at infrared wavelengths. Far- and near-ultraviolet 
fiuxes were obt ained from t he G alaxy Evolution Ex- 
plorer {GALEX, iMartin et aP l2005h data archive. These 
data are uncorrected for extinction and were assigned 
10 to 30 per cent uncertainties (larger than their quoted 
errors) owing t o this fact. Optical BVRI ph otometry 
was taken from [Bergeron. Leggett fc Ruij ll200ll'). supple- 
mented with ?7-band measurements in iMcCook fc SiorJ 
(200^), while ugriz photometry were available from the 
Sloan Digital Sky Su rvey Data Release 7 (SDSS DR7 
lAbazaiian et al.l I2OO9I '). Near-infrared JHK fiuxes were 
taken from the weighted aver age of photometry from 
'Bergeron, Leggett & Ruiz (2001') and the Two Micron All- 
Sky Survey (2MASS, Skrutskic ct al. 2006). The optical and 
near-infrared fluxes shown in Fig. [T] are all weighted equally 
with assumed 5 per cent errors. 

A thorough photo metric and trigonometric parallax 
analysis of GD356 by iBergeron. Leggett fc Rui3 (|200ll ) 
yields Tofi = 7510 K, logjQ((?/cms~^) = 8.14 and a helium- 
rich atmosphere for the DAe white dwarf. The top panel 
of Fig. [T] shows our flrst attempt to fit the data, repre- 
sented by points, with a model photosphere wit h both a 
7500K, logio(g/cms"^) = 8.0 DA spectral model (|Koesteil 
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Figure 1. Ultraviolet though near-infrared SED of GD356. The 
upper panel displays all available photometric data (filled sym- 
bols) for the white dwarf plotted beside a 7 500 K hydrogen-rich 
model (dashed line) and a 7 500K blackbody (dotted line). Circles 
are GALEX, upward triangles are UBVRI, downward triangles 
are ugriz and squares are JHK photometry. The lower panel 
shows fluxes for a 7 500 K helium-rich model (open circles) at the 
optical and near-infrared bandpasses, together with an 8 000K 
blackbody model (dotted line). This reproduces the helium at- 
mosphere model fluxes rather well. 



120091 ) and an identical temperature blackbody model. Nei- 
ther of these models is sufficient to account for the en- 
tire spectral energy distribution (SED) simultaneously. In 
the absence of a full spectral model for helium atmosphere 
white dwarfs cooler than 10 000 K, in the lower panel we plot 
the UBVRI and ugriz model fluxes as open circles for a 
helium-rich white dwarf of 7 500K and loRmf.q/cms"'^) = 8 
jHolberg fc Bere:eronll200^ : iFontaine. Brassard. Sz BergeronI 
I2001I ). against the photometry, represented as error bars. 
On top of this is plotted an 8 000 K blackbody model which 
mimics the 7 500 K helium-rich model fluxes quite well. We 
use this blackbody fit to extrapolate towards longer wave- 
lengths and iterate with the measured mid-infrared fluxes. 

Lastly we have analyzed archival Spitzer 
IRAQ images of GD 356 foll owing the methods of 
iFarihi. Jura fc Zuckermaiil (|2009l ). It is worth remarking 
that these images were free of potential flux-contaminating 
sources within the r — 3'.' 6 photometric aperture and that 
the signal-to-noise was sufficiently high that the flux errors 



Figure 2. Full SED of GD356. The upper panel shows all pho- 
tometry (filled circles) , now including the Spitzer IRAC flux mea- 
surements and the 8 000 K blackbody model, normalized to the 
same level as in Fig. [T] The lower panel is a linear plot of the 
infrared data and the potential modest excess at 8 ^m. 



at all wavelengths are dominated by IRAC calibration 
error. These previously unpublished fluxes are listed in 
Table [1] and plotted in the upper panel of Fig. [2] together 
with the shorter wavelength photometry and our selected 
model. The entire photometric SED, except possibly the 
7.9 /xm flux which appears to be in excess at the 2. Oct level, 
is fltted well by the 8 000 K blackbody model. The stellar 
image at this longest wavelength is both highly symmetric 
and sufficiently bright to show a faint Airy ring in the Q'.'6 
pixel" ^ mosaic but we are cautious about the interpretation 
of an excess at this level without more data (see the lower 
panel of Fig. [2| 



We foUow lFarihi. Becklin fc Zuckermaiil l|2008l ) and use 
the measured photospheric flux at 4.5 /im to place an 
upper limit to the mass of a possible substellar com- 
panion to GD356. For a white dwarf of mass 0.67 M© 
and effe ctive temperature 7 500K, the coo ling age is 
1.6 Gyr (jBergeron. Saumon fc Wesemael|[l995h . The main- 
sequence progenitor of GD 356 should have had a mass 
Mms ~ 3.25 Mp) according to the init i al-to-final mass re- 
lation (IWifliams. Bohe fc Koestej|2009l : iKahrai et allbOOSl : 
iDobbie et al.ll2006l ) and hence an estimated total lifetime of 
2.1 Gyr. According to models, at the d = 21.1 pc trigono- 
metric parallax distance to the white dwarf, an unseen 3a 
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(15 per cent) excess at 4.5 /^m places a companion upper 
mass limit of Afp < 12 Mj f or an age of 2.1 Gy r (Baraffe 
2007, private communication: iBaraffe et al.ll200a ). 

We note there are no whit e dwarfs with Spitzer-only 
exces s es owing to companions (iFarihi. Jura &: Zuckerman 



20091 : iFarihi. Becklin fc ZuckermanI l2008l : iMullallv et al 
20071) . All known infrared excesses from substellar compan- 
ions re veal themselves in t he near-infrared by 2 /xm at the 



latest ISteele et all [20091; IFarihi. Burleigh fc Hoard! 
Burleigh et al.1 l2006l: I Farihi. Zuckerman fc BecklinI 



2008 



20051 : 



Becklin fc ZuckermanI 1 19881 ). These results give strict up- 
per limits to the mass of Jupiter-sized companions to white 
dwarfs of typically between 5 and 20 M,j from young, less 
than 1 Gyr, to intermediate, 2 — 5 Gyr, total ages. The lim- 
its at GD 356 are therefore commensurate with other 4.5 
excess searches. 

The stringent upper limit that we have deduced for the 
mass of a possible secondary places new constraints on a 
binary accretion model for GD 356. First suppose that the 
substellar type companion is gaseous with Mp ^ 12 Mj. We 
expect such a companion to be tidally locked with the or- 
bital period just as in cataclysmic variables. For a Roche 
lobe filling gaseous companion of mass Mp orbiting a white 
dwarf of mass Mwd with separation a, the equivalent spher- 
ical Roche lobe radius Rl2 is given by 



Rl2 _ 2_( Mp 
a ^ 3I I Mp-hMwd 



l|Paczviiskilll97lf ). For Mp < Mwd Kepler's third law gi' 
Mo \ W P 



Rl2 = 0.05 



/ Mp Y ( py 



(1) 



(2) 



For 0.1 Mp/M j 70, the radi us of a planetary secondary 
is R2 ~ 0.1 -R0 (|Hubbardl I1994I ) . If such a secondary is to 
just fill its Roche lobe or lie within it the orbital period must 
satisfy 



(3) 



The limit of Mp < 12 A/j could therefore be consistent with 
a gaseous planet in a binary system with an orbital period 
PoTb ^ 2. 7 h in which the white d warf rotates at a period of 
115 min (|Brinkworth et al.l [2003 ) asynchronously with the 
orbit. Such a system could in principle be a post period 
bounce cataclysmic variable in which the companion mass 
has been reduced to low values by mass transfer. In such a 
case the white dwarf need not be tidally locked and should 
be rotating faster than the orbit because it is still accret- 
ing angular momentum, albeit at a slow rate. However, the 
time scale to reach post bounce perio ds of greater than 2.7 h 
would typically exceed a Hubble time (jKolb fc Baraffelll999l ) 
so that this possibility can actually be eliminated given the 
youth of the white dwarf. Alternatively a gaseous planet 
may have been dragged in to a close orbit during AGB evo- 
lution but not evaporated or it could have formed in a disc 
following a merger of two stars. 

However, for gaseous planets with Porb ^ 2.7 h, the 
mass transfer rate by Roche lobe overflow would be less 
than lO~^'^M0yr~^ when the orbital evolution is driven 
by gravitational radiation. For this white dwarf the corre- 
sponding specific accretion rate per unit area would be about 



10~^ gcm~^ s~^ if accreting material were to flow on to the 
10 per cent of the surface area of the white dwarf where 
a temperature inversion is seen. At these accretion rates a 
shock would not form and the atmosphere would be heated 
by particle bombardment and cooled by cyclotron emission. 
Deep heating would not be expected in the bombardment 
regime and heating to optical depths of order unity, as is 
required to explain t he emission lines in GD 356, is excluded 
(|Ferrario et al.|[l997h . Accretion from a companion at higher 
rates, perhaps by an irradiation induced wind, can also be 
rejected because the usual indicators of disc or funnel accre- 
tion are not seen in this system either. 

Eliminating these possibilities, we are once again 
led to favour the hypothesis of a conducting planet 
orbiting the magnetic white dwarf, with the excess 
emission po wered by the unip olar inductor mechanism 
(Li. Ferrario fc Wickramasinghd [1998). The planet must 
have a conducting core and be free of any atmosphere so 
that the inducted current dissipates in and heats the white 
dwarf atmosphere. At the mass of the Earth even a rocky 
planet in a close orbit would assume the shape of a Roche 
potential because tidal forces would be sufficient to melt 
any solid crust. For such a planet, of mean density of pp, 
the Roche limit obtained with equation ^ is 



> 47r 



1 



(4) 



This allows orbital periods of greater than about 4.7 h for 
pp = 5gcm~^. We have estimated the contribution that 
this planet would make to the observed energy distribution. 
A rocky planet with an albedo e in orbit at a distance a 
from the white dwarf, the equilibrium effective temperature 
is given, to a first approximation, by 



= (1 - e 



n1/4 / Rwd\^^^ , 

' \ 2a ) 



(5) 



For a bond albedo of 0.3, we find a planet temperature of 
560 K for an orbital period of 4 hours. Such a planet would 
contribute 1 per cent of the white dwarf flux at 7.9 pia if it 
had a radius of 1.4_R0. The 2a excess observed at 7.9 /xm 
could in principle be consistent with the presence of the 
hypothesized rocky planet, although other interpretations 
of this excess are also possible. 



3 DISCUSSION 

We must now ask how such a planet-like object can find itself 
so close to a white dwarf. We first show that it could not 
have been dragged in from a larger orbit by the progenitor 
of the white dwarf and then discuss the likelihood that it 
formed when a less massive companion white dwarf merged 
with GD 356. 



3.1 Planets orbiting white dwarfs 

The discoveries of Jovian planets orbiting evolved giant stars 
at distances of the order of 1 an demonstrate that such plan- 
ets can survive at leas t the early giant phases of evolution. 
iLovis fc Mma^ (|2007D estimate that at least 3 per cent of 
evolved giant stars (M ^ 1.8 A^q) have companions, in- 
cluding brown dwarfs, with A/pSini > 5Af,j. These results 



Does GD356 have a Terrestrial Planetary Companion? 5 



show that planets form in intermediate stars that evolve 
into white dwarfs. Some evidence that Jovian mass plan- 
ets may survive through the RGB and AGB phases to the 
white dwarf phase comes from the timing of pulsations in 
ZZ Get stars. It has been estimated that GD 66 may have a 
planet with Mp ^ 2.11 Mj orbiting at 2.4 au, although this 
is based on only one measured turni ng point of the orbi t 
and is therefore not well constrained (|Mullallv et al.ll2008h . 
Likewise V391 Pegasi, which is an extreme horizontal branch 
star, appears to have a planet wit h Mp ^ 3.2 Mj orbiting at 
about 1.7 au (|Silvotti et al.|[2007^ . However most SdB stars 
are in binaries so this particular system is likely to have been 
the end product of binary evolution. 

Gaseous planets that are initially close enough to in- 
teract with the expanding RGB or AGB star could simply 
be dragged in to a closer orbit by bow-shock and tidal drag 
during an ensuing common envelope phase of evolution or 
be completely destroyed by evaporation depending on their 
initial mass. The critical mass below which a Jovian planet 
is expected to evaporate in the envelope of the giant star 
before the envelope itself can be ejected is estimated to be 
about about ISMjup for a 1 Mq star but there are large 
uncertainties in this estimate related to parameters such as 
efficiency of common env elope ejection (jNelemans fc TaurisI 
ll998l : ISiess fc Livid[l999l ') . This estimate increases to 120Mj 
for a 5 Mq star. If the timescale for the common envelope 
phase were short enough it might be envisaged that a larger 
mass planet might have been partially evaporated down to 
about 10 Mj. However, given that evaporation should ac- 
celerate as the planet loses mass it is very unlikely that 
such fine tuning could have occurred. A main-sequence star- 
planet system that evolves through these phases would be 
seen either as a single white dwarf or as a white dwarf with 
a close planetary companion with a mass above this crit- 
ical value. GD 1400 and WD 0137 are close binaries with 
orbital periods of 10 and 2h re spectively ( (iMaxted et al.l 
l2006l : iFarihi fc Christop"hei1l2004l )). GD 1400 is a GO white 
dwarf, which is the remnant of AGB phase of evolution, 
while WD 1037 is a He white dwarf which has only passed 
through the RGB phase of evolution. The companions are 
50 — 60 Mj brown dwarfs. Thus the empirical evidence ap- 
pears to be that a companion of this mass can survive 
evaporation and in-spiralling during RGB/ AGB evolution, 
eject the envelope and be seen as a close binary now. We 
note, however, that I Villaver fc Livid jl999. ) have questioned 
whether such a close companion could survive the intense 
radiation to which it would be subjected to by the newly 
formed white dwarf and have proposed that such systems 
are more likely to have arisen from a merger of two white 
dwarfs. Against this hypothesis is the observation that main- 
sequence star-brown dwarf pairs appear to occur with about 
the same frequency as white dwarf-brown dwarf pairs indi- 
cating that the latter can form without the need for strong 
binary interaction. 

The fate of rocky planets, particularly that of the Earth 
itself, during the late st ages of stellar evolution has been 
looked at in some detail. ISackmann. Boothrovd fc Kraemen 
(Il993t ) thought that the Earth would escape because mass 
loss increases the orbital separati on faster than th e Sun 
grows on both the RGB and AGB. iRasio et al.l (|l99d ) went 
on to point out that tides, if a little stronger than ex- 
pected, induced by the Earth on the Sun might actually 



cause the Earth to spiral in during the RGB phase. Like 
ISackmann. Boothrovd fc Kraemen (1993) they worked with 
Reimers' ('Reimerd'l975') mass-loss rates and so claimed that 
survival of the RGB would lead to an orbit w ide enough to 
survive the AGB too. lRvbicki fc Denid (|200ll ') used more re- 
alistic thermally pulsing AGB models with weak mass loss 
on the RGB and so their stars experience many more ther- 
mal pulses and grow fast enough to swallow the Ea r th be - 
fore the end of the AGB. ISchroder fc Conan SmithI (120081 ) 
invoke even stronger mass loss towards the luminous tip of 
the RGB so that their Sun grows large enough to engulf 
a tidally spiralling Earth on the RGB but not the AGB. 
The orbital angular momentum of a low-mass terrestrial 
planet that is engulfed during the RGB phase is insuffi- 
cient to eject the common envelope. The orbit rapidly decays 
owing to tidal and bow-shock drag and the planet plunges 
into the central star and is destroyed. A planet that is first 
engulfed only during a thermal pulse on the AGB phase 
might have a better chance of survival. iRvbicki fc Denig 
(|200ll ') suggest that such a planet would be dragged in by 
10 — 70 Rp, per the rmal pulse and is also likely to be de- 
stroyed. IwHleT^Wu (2005 ) went back to t h e olde r AGB 
models of Sackmann. Boothrovd fc Kraemerl (|l993h to ar- 
gue that a fraction of such systems may survive stellar evo- 
lution and end up as close companions to the white dwarfs. 
However these models had fewer thermal pulses and did 
not grow as much on the AGB as is now thought. Indeed 
IWiUes fc W?3 (|2005l ) relied on the st ar shrinking sufficiently 
over t he last few pulses. Our mode ls (jStancliffe. Tout fc Pold 
l2004l : IStanclifre fc Glebbee^l2008D show AGB stars growing 
rapidly at the end of their lives after the onset of a super- 
wind. In which case engulfment cannot be avoided if a planet 
is dragged in. While uncertainty remains in the evolution of 
AGB stars it is possible that planets within a narrow range 
of separations from their stars might actually be dragged in 
to closer orbits. However the fine tuning of the orbital decay 
that would be required from pulse to pulse makes it very un- 
likely that the orbit could be reduced as much as required 
for GD 356. 

Assuming that any planet that might be sufficiently 
dragged in would be engulfed by the star's envelope, we 
must then consider the implicit assumption that an Earth- 
like planet is likely to evaporate during such periods of en- 
gulfment and so would not survive. To totally evaporate a 
planet, of mass Mp and radius i?p, it must absorb sufficient 
thermal energy to overcome its gravitational binding energy 



(,GMl 

Rr, 



10 



38. / Mp_ 



Rr, 



erg, 



(6) 



where ( — 3/5 for a uniform density body. For a typical 
rocky planet this is much larger than the total energy, of 
less than 10^^ erg, to sublimate, dissociate and ionize the 
planet. As long as the planet's core remains at a tempera- 
ture significantly lower than the ambient temperature bath 
the timescale for evaporation is just the Kelvin-Helmholtz 
timescale for the planet if it were radiating as a black body 
with the temperature of the bath. Thus the evaporation 
timescale 



(GM^/Rp 



(7) 



where T, is the temperature of the ambient giant stellar 
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Figure 3. The evolution of radius with time during the fifteenth 
thermal pulse of an initially 3 Mq star when it has a core mass of 
0.65 Mq. The star is beginning to lose mass in a superwind and so 
is expanding rapidly from pulse to pulse. The dashed line marks 
the maximum radius reached during the previous interpulse pe- 
riod. A planet at this radius is engulfed for more than 300 yr and 
reaches a depth of more than 50 Rq in the AGB envelope. 




Figure 4. The temperature structure with radius for the fifteenth 
pulse of our initially 3 Mq at t = in Fig. [3] when the star reaches 
its maximum radius. The dashed line marks the maximum radius 
reached during the previous interpulse period and so corresponds 
to the depth of a planet that just survived the fourteenth pulse. 



material. In equilibrium the thermal energy of the planet 
would exceed the gravitational binding energy when 

^^rcqMp>i5g„ (8) 

where R is the gas constant. The mean molecular weight 
/i depends on the composition and ionization state of the 
planetary material. At 10* K most constituents of the Earth 
are at least singly ionized and for a typical composition this 
gives fi ~ 15. At higher temperatures fi approaches two. 
Thus 

n^< 30000^^K. (9) 

Thus for an Earth-like planet at a depth that takes it to T, > 
30 000 K we may apply equation and tcvap is only about 
two days so the planet could not survive such conditions. 

Fig. [3] illustrates the evolution of the stellar radius of 
an initially 3Mq star with a current core mass of 0.65 
through a s ingle AGB pulse calculated w ith the Cambridge 
Stars code llStancliffe. Tout fc Pol3 120041 '). A planet that is 
located just outside the star prior to a pulse is engulfed by 
the expanding envelope for a period of more than 300 yr. 
Even if the planet were not dragged in, Fig. |4] demon- 
strates that it would be exposed to ambient temperatures of 
10,000 K for long enough to evaporate the more volatile el- 
ements. As the pulses proceed the mass of the planet would 
fall and the temperature and depth in the stellar envelope 
required for complete evaporation would rapidly diminish. 
We would not expect a rocky Earth-like planet to survive 
more than the first few pulses that engulf it. If the planet is 
dragged in it is exposed to yet higher temperatures for even 
longer. The temperature at the base of an AGB star's en- 



velope is in excess of 10^ K so a planet dragged into a close 
orbit within the giant envelope could not survive. 

For the extreme case of a planet engulfed at the very 
end of the thermally pulsing AGB when only about 0.05 Mq 
of stellar envelope remains the planet might survive partial 
immersion. However the binding energy of such an envelope 
is only about one thousandth of that required to bring the 
planet in to a close enough orbit around the white dwarf. 
We therefore conclude that a rocky. Earth-like planet cannot 
both survive evaporation and end up in a very close orbit. 

3.2 Merging double white dwarfs 

We consider a binary evolution scenario which leads to the 
formation of a CO white dwarf with a lower mass He or 
CO white dwarf companion through common envelope evo- 
lution. The two stars are subsequently drawn together by 
gravitational radiation and merge. During the merging this 
companion breaks up and forms a massive disc around the 
remaining CO white dwarf. Such a disc could be composed 
of CO-rich or He-rich material. The mass of the disc is 
somewhat more than one tenth of that of the accreting 
star and so becomes unstable to its own self gravity. The 
disc expands and cools as the central star accretes mat- 
ter. When the temperature in the outer disc is cool enough 
dust and rocks form and these clump to form a rocky plan- 
etary core. This model is very similar to that initially pro- 
pos ed to explain planetary co mpanions to millisecond pul- 
sars |Podsiadlgwskij|^ringle^^l^^ ( 1991. ). 

In order for the two white dwarfs to merge the common 
envelope process must leave them close enough for grav- 
itational radiation to act quickly. I Tout et all l|2008D have 
demonstrated how high magnetic fields in white dwarfs are 
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almost certainly generated in common envelopes from which 
the cores emerge already close together. Thus the high field 
in this system is evidence that the white dwarf most likely 
emerged from common envelope evolution with a close com- 
panion, in this case a second, less massive white dwarf. The 
mass of GD356 of 0.67 is already above the average 
for CO white dwarfs. In order to leave two white dwarfs 
that can merge the system must originally have been close 
enough that the evolution of both stars was curtailed by 
mass transfer. Thus we might envisage mass transfer from 
the initially more massive star to begin when it has a CO 
core of say 0.4 Mq . If this were to lead to a mild common 
envelope phase the orbit would then shrink so that the sec- 
ond star fills its Roche lobe early on red giant branch with a 
He core of about 0.3 Af0. Alternatively the first star might 
have filled its Roche lobe as a subgiant, evolved through an 
Algol phase to a helium white dwarf. Its rejuvenated com- 
panion could then go on to fill its own Roche lobe on the 
AGB followed by common envelop evolution that leaves its 
CO core in a close orbit. In either case for the final common 
envelope must leave the two cores sufficiently close that they 
can be driven together by gravitational radiation and their 
mass ratio must be less than 0.628 so that the ensuing mass 
transfer is dynamically unstable. 

When two white dwarfs merge, the more massive, being 
larger in radius, fills it's Roche lobe first. If the masses are 
sufficiently different stable mass transfer could follow with 
the orbit widening as the mass-losing star grows in radius. 
However if the donor is more than 0.628 times the mass 
of the accretor the mass transfer is unstable because the 
white dwarf grows faster than its Roche lobe expands. In this 
case numerical simulations show that the less massive white 
dwarf is indeed tidally disrupted and accret ed on to the more 
massive in through a t hick accretion disk CilBenz et al ] |l990l : 
[Cucrrcro J., Garci'a-B erro E. fc Isern J.I [20041 '). So that the 
natural outcome is a hot white dwarf surrounded by a thin 
remnant disc that contains most of the angular momentum. 
The nature of this disc is likely to be very unusual, given that 
white dwarfs are composed primarily of carbon and oxygen, 
rather than the hydrogen and helium of more traditional 
circumstellar discs. 

The formation of planets in such discs around neutron 
stars and white dwarfs has been discussed bv lHansenI (|2002l ') 
and iLivio, Pringle fc WoodI l|2005l ). The expected outcome 
depends rather critically on the viscosity of the disc. Ini- 
tially it has an outer radius of 10^ — 10^^ cm, determined by 
the orbital angular momentum of the disrupted companion. 
As accretion proceeds the disc expands and cools. While the 
viscosity is determined by the magneto-rotational instabil- 
ity it is strong only when the disc is ionized and the viscos- 
ity is negligib le outside this region. If sufficient gas persists 
iHansenl (|2002l ') then finds that planets form in the quiescent 
outer disc but at a higher temperature than for a hydrogen 
rich composition, because of the higher ionization potential 
of carbon and oxygen, by processes similar to those that are 
believed to have occurred in the early solar system (|Lissaueij 
1 19931 ). They predict that planets of 30 - 300A/e form in CO 
rich discs and are located within 0.2 au. A similar scenario 
is appropriate for He rich discs. Planet formation may take 
10* yr. The temperature of the white dwarf indicates that 
it has been cooling in excess of 10^ yr which easily accom- 
modates this along with sufficient time for any remnant disc 



to disperse. We might further speculate that all volatiles, 
and perhaps some transitional elements, such as magnesium 
and silicon, may evaporate and be lost in the very near en- 
vironment of a hot and relatively luminous, newly merged 
WD, while refractories, such as calcium, titanium and alu- 
minium, would more likely be retained, perhaps leading to 
a refractory-metal planet. In particular much of the oxygen 
that would otherwise form oxides would almost certainly 
evaporate and be blown away by radiation pressure so that 
metals that would otherwise oxidize would be left to form a 
substantial metallic core. Such a planet might be more like 
Mercury than the Earth in composition. 

Such a planet would have a metallic core of mass 0.03 — 
IMe (for a CO rich composition). Once within 10 Rq of the 
highly magnetic white dwa rf orbital energy can drive the 
unipolar induction current fci. Ferrario fc Wickramasinghd 
1998). Owing to its atmosphere the planet's effective resis- 
tivity is initially larger than that of the white dwarf so that 
energy, which is extracted from the orbit, is mainly dissi- 
pated in the planet during the early phases of the mag- 
netic interaction. As the planet drifts in, this heating facil- 
itates the evaporation of its atmosphere, until the effective 
resistivity of the planet becomes smaller than that of the 
white dwarf's atmosphere. The heating then occurs mainly 
in the white dwarf atmosphere and th e model presented by 
iLi. Ferrario fc Wickramasinghd (|l998l ) for GD 356 becomes 
applicable. 

Some observational support for the possibility of the 
formation of a second generation of planets around a star 
that is evolving into a white dwarf has been provided by 
the discovery of a substantial classical T Tauri type dust 
disc surrounding an accreting firs t giant ascent giant star 
TYC 4144-329-2 (|Melis et al.ll2009l 'l. It has been speculated 
that the observed disc may have resulted from the common 
envelope interaction with a low mass stellar or substellar 
companion. 



4 CONCLUSIONS 

We have presented an analysis of archival Spitzer observa- 
tions of GD 356 that shows that a firm upper limit of 12Mj 
can be placed on the mass of a possible companion. The 
new observations places further constraints on the orbital 
parameters of a possible binary companion. In agreement 
with previous investigators, we have again argued that ac- 
cretion heating due to mass transfer from a companions is 
an implausible mechanism for explaining the anomalous line 
emission in this star. The unipolar inductor model, which re- 
quires GD 356 to have a rocky planet with a metallic core as 
a companion, remains the best explanation for the peculiar 
properties of this unique star. 

Theoretical estimates indicate that it is unlikely that an 
Earth type planet that was orbiting the main-sequence pro- 
genitor of a white dwarf, at a distance that would allow it to 
be engulfed by the expanding envelope during RGB/ AGB 
phases of evolution, would survive the subsequent evolution 
of the parent star and be seen as a close companion to the 
white dwarf. Such a planet would either be evaporated as it 
is dragged into the core of the star during the RGB/ AGB 
phases of evolution or be left in an orbit at a much larger 
radius of several hundred solar radii. If the unipolar indue- 
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tor model for GD 356 is to be viable an alternative origin 
must be sort for its close companion. We have argued that 
the planet probably formed from in an accretion disc follow- 
ing the disruption of a white dwarf in a rare merger event. 
GD 356 would thus be the white dwarf counterpart of the 
millisecond binary pulsar PSR1257+12 which is known to 
host a planetary system. 

In conclusion, we note that the unipolar inductor model 
for GD 356 makes a definite prediction that could be veri- 
fied by future observations. Two fundamental and distinct 
periods, the rotation period of the white dwarf and the or- 
bital period of the planet around the white dwarf, are ex- 
pected to be seen in the emission line fiux or in any com- 
ponent of the continuum flux attributable to the heating. 
iLivio. Pringle fc Safferl l|l992f ) argue that planets are most 
likely to be found around massive white dwarfs because these 
are more likely to have merged in the past. We would add 
that the white dwarf should not only be massive but also 
posses a high magnetic field, created during the common 
envelope evolution that must have preceded the merging. 
It is around such white dwarfs that the search for planets 
should be concentrated and, if metallic, such systems might 
also show up as unipolar inductors. 
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